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We extensively investigate two-step shape invariance in the framework of AA-fold supersymmetry. 
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1 -Q . general form of type A 2-fold supersymmetry, we systematically construct two-step shape-invariant 
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I. INTRODUCTION 

In the formalism of supersymmetric quantum mechanics (SUSY QM) [1, 2] potentials 
with unbroken SUSY and shape invariance (SI) [3] can be exactly solved by a well-known 
standard procedure [4, 5]. The potential algebras of these systems have also been identified 
[6-9] providing an alternate method of finding exact solutions. Mainly three classes of SI have 
been studied in the literature, namely, (1) translational class [10-13] - when the parameters 
differ only by a constant i.e. b^^' = fe^*^-* + a, a being a constant, (2) scaling class [14, 15] - 
when the parameters are related by a scaling factor i.e. 6*-^^ = qb^^^ with < g < 1, and (3) 
exotic class [14] - when the parameters are related by other equalities such as b^^^ = q(b^'^^y 
and b^^^ = qb^^y{l +^6'-°^). It should be mentioned here that SI is neither a necessary nor a 
sufficient condition for exact solvability. In fact, some exactly solvable potentials are shown 
not to be SI [15] . Recently a complete set of additive SI potentials have been generated from 
an Euler equation [16]. The SI condition is also studied in the context of fractional SUSY 
QM of order k (/c = 3, 4, . . .) generalizing the Z2 grading of the relevant Hilbert space to a 
Zk grading [17-20]. 

The idea of SI can be extended to the more general concept of SI in two and even multi- 
steps. While there are quite a number of works on one-step (namely, ordinary) SI, the 
literature dealing with SI in two or more steps is rather scanty. In [14] the authors first in- 
troduced the concept of SI in two steps to enlarge the class of exactly solvable Hamiltonians. 
The two-step SI approach was utilized for dealing with SUSY QM problems with sponta- 
neously broken SUSY [21]. For a quantum mechanical system with position-dependent mass 
the same approach was used to handle broken SUSY problem [22]. Recently a class of solv- 
able potentials of translational SI in two steps were obtained [23, 24]. It was found that 
discontinuity at some points was a characteristic of the two-step superpotentials, therefore 
giving rise to Dirac delta-function singularities in the corresponding potentials if they are 
considered in the whole line x G (— C)0, +00). The translational SI potentials were shown to 
possess a potential algebra involving three generators of angular momentum type. The po- 
tential algebra for the case of SI in k steps {k being an arbitrary positive integer) was found 
[25] to be equivalent to the generalized deformed oscillator algebra that had a built-in Z4 
grading structure and was constructed in terms of the generators of the deformed harmonic 
oscillator (J, A, A\ N) as well as the grading generator T of the cyclic group of order k. The 
obtained potentials included the cyclic SI potentials of period A; as a special case. 

On the other hand, the framework of A/'-fold SUSY [26-28] has gained much interest 
during the last few years. It is possibly one of the most powerful frameworks among various 
existing methods for finding or constructing quantum systems which admit exact solutions 
in a certain sense. In particular, a few remarkable features are in order here: 

1) A/'-fold SUSY includes all the ordinary SUSY methods as its special cases and is equivalent 
to weak quasi-solvability [27] which is less restrictive concept than quasi-exact solvability 
[29, 30]. 

2) An arbitrary one-body quantum Hamiltonian which admits two (local) solutions in closed 
form belongs to a special class of A/'-fold SUSY, namely, type A 2-fold SUSY [31, 32] irre- 
spective of whether or not it is Hermitian, VT symmetric [33], pseudo-Hermitian [34] and 
so on. 

3) Many of the so far constructed (quasi-)solvable position-dependent mass quantum Hamil- 
tonians are also realized as special cases of type A A/'-fold SUSY [35] . 

We note that almost all the models having essentially the same symmetry as A/'-fold SUSY 



but called with other terminologies in the literature, such as Poschl-Teller and Lame po- 
tentials, are actually also particular cases of type A A/'-fold SUSY. For a review of A/'-fold 
SUSY, see Ref. [36]. Recently in Refs. [37, 38], the necessary and sufficient conditions for 
type A A^-fold SUSY Hamiltonians to possess intermediate Hamiltonians were derived for 
M = 2 and 3. As a by-product, some well-known translational SI potentials were naturally 
obtained. It indicates that one can in principle go one step further and consider two-step SI 
also in the framework of A^-fold SUSY. 

In this paper, we extensively investigate two-step SI in the framework of A/'-fold SUSY. 
A key ingredient of our approach relies on the observation that any two-step SI system pos- 
sesses type A 2-fold SUSY with an intermediate Hamiltonian whose analytic and algebraic 
structures are both well understood. The latter fact indeed enables us to make a sophisti- 
cated analysis on general aspects of two-step SI without recourse to any specific assumption 
or model. Employing the general form of type A 2-fold SUSY, we systematically construct 
two-step SI potentials without relying on any ad hoc ansatz. In addition to the well-known 
ordinary SI potentials, we successfully obtain several novel two-step SI ones which have not 
been reported in the existing literature. Furthermore, some of the latter two-step SI are 
conditional and thus the corresponding systems are conditionally solvable. 

We organize this paper as follows. In Section II, we review the concept of two-step SI 
and discuss its general aspects. In particular, we show that any two-step SI system has 
type A 2-fold SUSY with an intermediate Hamiltonian and second-order paraSUSY and 
that two-step SI always means solvability and actually lies between solvability and ordinary 
SI. In Section III, we investigate two-step SI for particular models which can be realized 
in the cases of type A A/'-fold SUSY for Af > 2. The obtained two-step SI potentials 
include the well-known ones such as (radial) harmonic oscillators, Morse, Scarf, Poschl- 
Teller potentials, and so on. In addition, we also obtain some novel potentials which are 
expressible in terms of elliptic, exponential (including hyperbolic and trigonometric), and 
rational functions. In Section IV, we construct more general two-step SI potentials under a 
less restrictive condition. We obtain more novel two-step SI potentials some of which admit 
analytic expressions only implicitly. Most of the new potentials obtained in Sections III 
and IV do not have ordinary SI and some of their two-step SI are conditional. Finally in 
Section V, we summarize the results and discuss perspectives and future issues. 

II. TWO-STEP SHAPE INVARIANCE AND TYPE A 2-FOLD SUSY 

First of all, let us introduce the concept of two-step SI. Suppose we have a system com- 
posed of two sets of SUSY QM 

2Hf{c) = Af{c)Af{c), Af{c) = T^ + W,{x;c), (z = 0,l), (2.1) 

where c denotes a set of all the parameters involved in the system, which satisfies the 
following constraint 

H^{c) = H^ic) + R{c), (2.2) 

where R{c) is a constant depending only on c. According to Ref. [14], the system is said to 
have two-step shape invariance if Hq and H^ satisfy 

Ht{c^'^)=H^{c^'^) + RM'^), (2.3) 



where i?2(c^°^) is another constant and c*^^^ = c'^^^(c*^°^) is another set of parameters both of 
which depend only on c^'^'. 

To begin with, we shall show in what follows that any such a system has type A 2-fold 
SUSY with an intermediate Hamiltonian studied in Ref. [37]. We first note from (2.1) and 
(2.2) that Hq{c) and H^{c) + R{c) are intertwined by the second-order linear differential 
operator A^^Ag^ as 

A-,{c)A^{c)H^{c) = {Ht{c) + R{c))A-,{c)A-,{c), (2.4) 

and thus form a pair of 2-fold SUSY. It is evident from (2.2) that the latter 2-fold SUSY 
system possesses an intermediate Hamiltonian. In addition, we can have an expression of 
the kernel of Ai{c)Aq{c) in closed form as 

kerA^ (c)A^(c) = i^, I dx exp( f dx" {Wq{x'- c) - Wi{x"; c)) 

xexpf-/ dx'Wo{x';c)\ (2.5) 



which is preserved by Hq{c) due to the intertwining relation (2.4). Hence, Hq[c) admits 
two linearly independent local solutions in closed form. The latter fact is a necessary and 
sufficient condition for a one-dimensional Hamiltonian to belong to type A 2-fold SUSY 
proved in Ref. [32]. Therefore, we finally conclude that any two-step SI system has type A 
2-fold SUSY with an intermediate Hamiltonian. 

It is now clear from the above process of the proof that any two-step SI system (2.1)- 
(2.3) has one-to-one correspondence with a type A 2-fold SUSY system with an intermediate 
Hamiltonian {H^,H'\P^~ = P21P22) in Ref. [37] by 

H^ ^ H-, H+ ^ H'\ H+ + R^ H+, A^ o P2I, M ^ P^- (2-6) 

To investigate general structure of two-step SI systems in the framework of type A 2-fold 
SUSY with an intermediate Hamiltonian, it is convenient to make a 'gauge' transformation 
from the physical x-space to a gauged 2;-space defined by 

H^ = e^^H^e-"^^, P^ = e^^ P:^^-'^^ ■ (2-7) 

According to Ref. [37], the general form of a gauged type A 2-fold SUSY system with an 
intermediate Hamiltonian is given by 

H-ic) = -A(:)^ - Q{:\ 6)^ + \q'{:; b) - R, (2.8) 

H*{c) = -A(z)^ - Q{z; b)± - |q'(.; b) + ^<^i||p> - R, (2.9) 

A-W = (.T^. P+(c)=e-/-4AV-5 = (/)^(iL+«i|5)y. (2.10) 

with the set of parameters c = {b,R}, where A(z) is an arbitrary function of z at present 
and Q{z; b) is a polynomial of at most first-degree in z 

Q{z;b) = biz + bo. (2.11) 



The function z' = dz/dx which connects the variable z in the gauged space and the coordi- 
nate X in the physical space is determined by 

^'(^)' = 2^WL=.(.)- (2-12) 

It can be easily checked by a direct calculation that if ^ are intertwined by P^ as 

P^{c)H-{c) = H^{c)P^{c), P^{c)H^{c) = H~{c)P+{c). (2.13) 

With a factorized form of the 2-fold supercharges 

P^{C) = P2"l(c)P22(c), P^C) = Pt2{c)P^l{c). (2.14) 

the three component Hamiltonians H~ , W^ , and H^ are written as 

2H~{c) = P}^{c)P^^{c) + 2C22(c), (2.15) 

2H'\c) = P^^{c)P+^{c) + 2C22(c) = P2^(c)P2l(c) + 2C2i(c), (2.16) 

2#+(c) = P2i(c)P2"^(c) + 2C2i(c). (2.17) 

The type A gauged 2-fold supercharges (2.10) admit a factorization (2.14) with 



f. , ^ / / d A'iz) \ ~ , , , d 

P22(C)- ^^^^+ 2A(^) J, ^llCJ- ^^^^+ ^(^) 



(2.18) 



It is easy to check that the pair of Hamiltonians H^ (2.15) and (2.17) obtained from (2.18) 
actually coincides with the pair of type A gauged Hamiltonians (2.8) and (2.9) with 

2C22 = &i - 2R, 2C21 = -hi - 2R. (2.19) 

The intermediate gauged Hamiltonian H^^ (2.16) with respect to the factorization (2.18) is 
calculated as 

As was shown in Ref. [37], the type A gauged 2-fold supercharges (2.10) admit, in addition 
to (2.18), a one-parameter family of factorizations, denoted by symbols with wide hat, with 



dz 2A{z) z + zqJ' \dz z + zq 

^''~ ' [dz^ 2A{z) ^ z + zoj' ^ ^ 



P+ - -z' 

-' 21 ~ -^ 



d , Q{z;b) 1 



dz A(z) z + zq 



where ^o G C is a parameter. Then, a type A 2-fold SUSY system admits another inter- 
mediate gauged Hamiltonian iP^ with respect to one of the above one-parameter family of 



factorizations if and only if 61 7^ 0. In fact, if it is the case, the pair of gauged Hamiltonians 
H^ (2.15) and (2.17) calculated with the second factorization (2.21) coincides with the pair 
of type A gauged Hamiltonians (2.8) and (2.9) with 

2(522 = -hi - 2R, 2C21 = 61 - 2R, (2.22) 

provided that Zq = 60/^1 • The second intermediate gauged Hamiltonian if' ^ with respect to 
the second factorization (2.21) defined via (2.16) with the corresponding hatted quantities 
is calculated as 

Another remarkable aspect of type A 2- fold SUSY systems with an intermediate Hamiltonian 
is that they possess second-order paraSUSY as well [37]. Introducing a triple of operators 
(i/p, Qp) as 

if P = H-{i,^f{i;^f + W\i^ti^^ - {i^t)\^^?) + H\i;m^^)\ .2 24) 

Qp=^2+2(^?)'^P+^2^V'It(^p)^ Q+ = P2"2^p(^+)2 + P2'l(V^lt)'^F, 

where H^ are a pair of type A 2-fold SUSY Hamiltonians, W^ is (one of) its intermediate 
Hamiltonians, and il)-p are second-order parafermions satisfying 

(^^)V0, (V^^)' = 0, {^p,V'it} + {(V'p)',(V^i^)'}=2/, (2.25) 

we see that the triple (ifp, Qp) satisfies the second-order paraSUSY relations in [39] 

(Q^)Vo, (Q^)' = o, [Q^,ifp]=o, 
(Qp)'Q? + QtQ^Qt + Q?(Qp)' = 4Q^^p, 

as well as the generalized 2-fold superalgebra 

{Q^fiQt? + Qp(Q?)'Qp + {QtfiQ^f = KH^f - {h)\ (2.27) 

where hi is the same parameter as the one appeared in (2.11). 

To see what are significant consequences of two-step SI, let us first show that if a pair of 
type A 2-fold SUSY gauged Hamiltonians H^ satisfies the two-step SI condition (2.3), that 
is, 

H+{c^^^)) = H-{c^^^^^)) + i?2(c(2'=)), (2.28) 

where we have generalized c*^"^ and c^'^^ in (2.3) to c^^*'^ and c^^'^"'"^) [k = 0, 1, 2, . . .), and if, 
in addition, H~ {c^'^^^'^'>) preserves a linear space V~(c(^'^"'"^)), then H~{c^'^'^'>) preserves the 
space ig (c'^^'^^)V~(c'^^'^"^^)). In a mathematical language, we have 

^-(c(2fc+2))^-^^(2/c+2)) C V^(c(2'=+2)) 

=^ ^-(c(2^-))P2+(c(2^))V^(c(2^+2)) ^ P2^(c(2'=))V-(c(2^-+2)). (2.29) 



Indeed, it follows from (2.14)-(2.17) that 

= P2+2(C^''^) te(c(2'))P2~l(c^''^) +2C2l(c(2^')) P2"i(c(2'=))V-(c(2'=+2)) 

= P2t(c^''^)^2"l(c^''^) fel(c^''^)^2'\(c^''^) + 2C2l(c(2'=')] V-(c(2'=+2)) 
= 2P2^(c(2*^))^+(c(2^))V-(c(2'=+2)). 

Hence, under the condition (2.28) we obtain 

i^-(c(2^))P2+(c(2'=))V-(c(2'=+2)) ^ p+(c(2^-)) [if-(c(2'=+2)) +P2(C(''^)1 V~ {c^^''^'''^) 



(2.30) 



and thus conclude (2.29). On the other hand, it is evident from (2.13) that H preserves 



the kernel of P2 which is at most two dimensional: 



H~{c)%{c) C %{c), %{c) = kerP2^(c) = (l,^. 



(2.31) 



and thus it is quasi- solvable [36]. Combining (2.29) and (2.31), we see that the Hamiltonian 
H~(c^'^^) preserves an infinite flag of finite-dimensional linear spaces 



V,^(c(o))c V7(c(°))c---CV, 



2n\ 



,(0)^ 



C 



(2.32) 



where V2„(c(°)) is defined by 



V. 



2n\ 



,(0)^ 



V2-(CW)+P2+(CW)V2-(C(^)) + --- 
■ ■ ■ + P,+ {C^'^)P,+ {C^^^) ■ ■ ■ P2+(c('""'^)V2^(c('""'^) 



(2.33) 



Hence, any two-step SI Hamiltonian is solvable [36] unless there is a natural number n such 
that V2n{c^^^) C V2^(c*^*^)). We also note that if an intermediate Hamiltonian W^ satisfies 



^il(c(2^'+l)) = iJ-(c(2A^+2)) ^ ^^(c(2fc+l))^ 



(2.34) 



where Ri(c^'^''~^^^) is a constant depending only on c^'^''+'^) ^ then the three Hamiltonians H 
Wil, and H^ constitute a sequence of ordinary SI with 



^+(c(2'=)) = ^il(c(2fc+l))+p^(c' 



(2fc)^ 



(2.35) 



where R2{c^'^''^) = Ri{c^'^''~^^^) + Pi(c'^^'^)). Conversely, it is evident from the definition that 
any ordinary SI Hamiltonian is two-step SI. It means that ordinary SI can be regarded as a 
sufficient condition for two-step SI. Summarizing the above arguments, we have the following 
relation: 



(Solvability) D (Two-step SI) D (Ordinary SI). 



We shall say that two-step SI is irreducible if it is not ordinary SI simultaneously, and is 
reducible otherwise. 

The pair of type A 2-fold SUSY potentials V^ with an intermediate Hamiltonian is given 

by [37], 



V^(x,c) 



A{z)A"{z)-^A'{z)^-{b^z + boy 



AA{z) 
T2{{b^z + bo)A'{z)-2b,A{z)) 



- R, (2.36) 

z=z{x) 



where the set of parameters is c = {bi, bo, R}. The potential terms in the first and second 
intermediate Hamiltonians H^^ and H^^, the latter of which exists only when 61 7^ 0, are 
calculated as 



ya..__ ^"W A'{zr-A{b^z + boy 



4 16A{z) 



- R, (2.37) 

z=z(x) 



ri2r^ „^ ,nu^ „^ hA'{z) , 2{b,)'A{z) 



r^(x;c) = r^(x;c)--^^ + 



biz + bo {biz + bo) 



- R. (2.38) 

z=z{x) 



In what follows, we put i? = to remove an irrelevant parameter without any loss of 
generality. From the form of a pair of potentials (2.36), we immediately see that there is 
always two-step SI of the following type 

V-^{x;bi,bo) = V-ix;-bu-bo), (2.39) 

that is, the relation (2.3) holds with bf^ = -bf\ bf = -6S,°\ and ^2(&f\C) = 0> 
irrespective of the form of A{z). The latter two-step SI usually changes the normalizability 
of the corresponding solvable sectors and thus can be employed to convert a solvable model 
without normalizable eigenfunctions into an exactly solvable one. A similar tactics was 
already demonstrated, although with a different type from (2.39), in Ref. [21]. We shall 
say that the two-step SI under consideration is reflective if it is characterized by (2.39). In 
general, reflective two-step SI is irreducible. Our central concern is now what kinds of A[z) 
can admit non-refiective two-step SI. 



III. POLYNOMIAL A{z) 

To begin with, we recall the fact that in type A A/'-fold SUSY for M > 2 the function 
A[z) must be a polynomial of at most fourth degree: 

A{z) = a^z^ + a-sz^ + a2z'^ + aiz + ao- (3.1) 

Hence, it is natural to examine first the case of this polynomial A{z) of at most fourth 
degree. In the followings, we present the forms of potentials resulting from a polynomial 
A{z) of from zeroth to fourth degree. 

Case 1-1: 04 = 03 = 02 = ai = & oq 7^ 

Let us first consider the case of a zeroth degree A{z). In this case, the pair of potentials 
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(2.36) are given by 



V^{x;h,bo) 



{hfz^ + 2biboz + {bofT^aobi 



4an 



z=z{x) 



In addition to (2.39), they have non-reflective two-step SI as 

V+ix;bi,bo) = V-ix;bi,bo) -2b,, 



(3.2) 



(3.3) 



that is, the relation (2.3) holds with bf^ = bf\ bf = 6[,°\ and R2{bf\b'^^) = -2bf . The 
first intermediate potential (2.37) is calculated as 



V'\x-M.h. 



)i)V + 26i6o^ + (M' 



4an 



(3.4) 



z=z(x) 



and thus it together with V^ constitutes a sequence of ordinary SI as 

v'\x-MM) = v-{x-MM-K 



(3.5) 



Therefore, the system has reducible two-step SI. On the other hand, the second intermediate 
potential (2.38) which exists only when 6i 7^ reads as 



v'\x-MM = v'\xMM + 



2ao(6i)2 



ipiz + b^y 



(3.6) 



z=z(x) 



and has no SI with the other potentials. The equation (2.12) in the present case is integrated 

as 



z = \/2an X. 



(3.7) 



Substituting (3.7) into (3.2), (3.4), and (3.6), we finally obtain the potentials in the x-space. 
The system is well-known harmonic oscillators. 

Case 1-2: 04 = 03 = 02 = & ai 7^ 

Next, we shall consider a first degree Ai^z^ with 04 = 03 = 02 = and a^ 7^ 0. In this 
case, the pair of potentials (2.36) are given by 



v'^ix-MM) 



ihr 



z + 



bittQ 26o±3ai\ fbittQ 26o ± oi 



4ai A{aiZ + ag) \ o-i 

- - — 26o ± 2ai 

4ai \ ai 



Oi 



z=z{x) 



In addition to (2.39), they have non-refiective two-step SI as 

V+{x- 61, 60) = V-{x; 61, 60 + 2ai) - 26i, 



(3.8) 



(3.9) 



(2) _ , (0) , (2) _ , (0) 



that is, the relation (2.3) holds with b\ = b\ , 6q 
The first intermediate potential (2.37) is calculated as 



(0 >(o)^ 



^0 



2ai, and R2{bi , 69 



-2b'; 



(0) 



nU^.r^ .^ (K 



V''ix;b^,bo) 



-z + 



61 ao 26o + ai\ fbiao 26o — cti 



4ai 4(ai2; + ao) \ '^1 

61 fbiao 



ai 



4ai V ai 



2br 



(3.10) 



2=z(x) 



and thus it together with V^ constitutes a sequence of ordinary SI as 

V'^{x] bi,bo) = V~{x; fei, 60 + ^i) - h, 
V+{x; bi,bo) = V'\x; bubo + ai) - b^. 



(3.11) 



Therefore, the system has reducible two-step SI. On the other hand, the second intermediate 
potential (2.38) which exists only when 61 7^ reads as 



v'\x-MX) = v'\xMM + 



aibi 



+ 



aibo)bi 



biz + bo {piz + 60)^ 



(3.12) 



z=z{x) 



and has no SI with the other potentials. The equation (2.12) in the present case is integrated 
as 



Z = — X . 

2 ai 



(3.13) 



Substituting (3.13) into (3.8), (3.10), and (3.12), we finally obtain the potentials in the 
x-space. The system consists of well-known radial harmonic oscillators. 

Case 1-3: 04 = as = & 02 7^ 

Next, we shall consider the case of a second degree A{z). In this case, the pair of potentials 
(2.36) are given by 



V^(x;bi,bo) 



(&i)' + (o2)' , (20260 -ai&i)(&i± 202) 



+ 



4a2 4:a2{a2Z^ + aiz + ao) 

02(260 ± ai)(26o ± 3ai) — 4ao(6i ± (3.2)(6i ± 3a2) 



z=z{x) 



1602(02^^ + aiz + ao) 
In addition to (2.39), they have non-reflective two-step SI as 

V^{x; 61, 60) = V-{x] 61 + 4a2, 60 + 2ai) - 2(6i + 2a2 

that is, the relation (2.3) holds with bf^ = bf^ + 4a2, 6?^ = 6^ + 2ai, and i?2(6r,&o 
—2(b[ + 202). The flrst intermediate potential (2.37) is calculated as 



(3.14) 

(3.15) 
(0 ^(o)^ _ 



V'\x-b,,k 



(61)^ + (02)^ (20260 - Qi6i)6i 

4a2 4:a2{a2z'^ + aiz + ao) 

02(260 + ai)(26o - fli) - 4ao(6i + a2)(6i - 02) 



1602(02^^ + o-iz + ao) 



(3.16) 



z=z{x) 
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and thus it together with V^ constitutes a sequence of ordinary SI as 

V'^{x; bi, bo) = V~{x; bi + 2a2, &o + cli) - bi - 02, 
V~^{x; 61, 60) = V'^{x; bi + 2a2, &o + ai) - &i - 02- 



(3.17) 



Therefore, the system has reducible two-step SI. On the other hand, the second intermediate 
potential (2.38) which exists only when 61 7^ reads as 



V''{x;b,,bo) = V'\x;b,,bo) + 



aibi - 2a2bo ao{bi)'^ - aibibp + a2{boY 
biz + 60 



{biz + bo) 



(3.18) 



z=z{x) 



and has no SI with the other potentials. The equation (2.12) in the present case is integrated 

as 



2a2 2 



Oo 



2a2 8(02)2 



/2a2 X 



(3.19) 



Substituting (3.19) into (3.14), (3.16), and (3.18), we finally obtain the potentials in 
the x-space. The system consists of well-known Morse, Scarf, or Poschl-Teller poten- 
tials, depending on the relative relation among given values of the parameters Oj {i = 0, 1, 2). 

Case 1-4: 04 = & 03 7^ 

Next, we shall consider the case of a third degree A{z). In this case, the pair of potentials 
(2.36) are given by 

.±/ , 3a3 02 ±861 3a3(3ai =F 8&0) - 0,2(302 =F 8&1) - 4(61 )2 ^ 



v^ix-MM 



16"^^ 16 lQA{z) 

27a^ao — ai(3a2 =F I6&1) — 860(^1 ± 202) 

lQA{z) 
ao(9a2 ± 246i) - ai(3ai ± 860) - 4(6o 



lQA{z) 
The first intermediate potential (2.37) is calculated as 



(3.20) 



z=z(x) 



v'\x-MM) 



15a3 5a2 4(6i)2 + 30301 - (02)^ ^ 
~l^^ ^ Te" ^ lQA{z) ^ 

86ifeo + 90300 — 02O1 4(60)^ + 302O0 — (ai)2 

H ,^ ,, , z H 



lQA{z) 



lQA{z) 



(3.21) 



z=z{x) 



On the other hand, the second intermediate potential (2.38) which exists only when 61 7^ 
reads as 



+ 2 



0360 oi(6i)2 - 202&160 + 303(^0 
61 bi{biz + bo) 

aoibi)^ - ai{bi)%o + 02^1 (&o)^ - 03(^0)^ 



61(612 + 60)2 



(3.22) 



z=z(x) 



Although V^ and V have the same functional dependence on the variable z, it does not 
necessarily mean that there is a transformation of the parameters 61 and 60 which lead to 



11 




non-reflective two-step SI. In fact, they satisfy the condition for two-step SI (2.3) if and only 
if 

-(6^ + 202) 603 

60 + 4ai b^ -4a2 \ { IX ] = 0, (3.23) 

6ao b^ — 2ai 

where bf = bi ± b\ . The trivial solution b^ =bQ = f) just produces reflective two-step SI 
(2.39). To have a non-trivial solution, the parameters must satisfy 

603(60 + 4ai) = -(&r + 2a2)(6r - 4a2), (3.24) 

(6r + 202) (60 - 2ai) = -36a3ao. (3.25) 

We flrst consider the particular case when b^ = — 2a2- In the latter case, the set of equations 
(3.24) and (3.25) reduces by the assumption 03 7^ to 

b+ = bQ + Aai = ao = 0. (3.26) 

Hence, we have 

A{z) = asz^ + a2Z^ + a^z, bf^ = bf^ - 2a2, &[,^^ = -&[,°^ = -2ai. (3.27) 

In this case, the pair of potentials (3.20) reads as 

,± ^ , _ 3a3 ^ , 02 ± 861 3a3(3ai t 860) - 02(802 T 861) - ^bif 



V [x; 61, bo) = —-Z + — — — — ■ ^ z 

16 16 Ibla^z-^ + a2Z + ai) 



ai{3a2TlQbi) + 8bo{bi±2a2) ai (3ai ± 860) + 4(6o) 



16(03^^ + a2Z + ai) 16(03^^ -|- a2Z + ai)z 

We can easily check that the latter pair satisfles 



(3.28) 



z=z(x) 



V^{x- 61, bo)A{z) = V-{x- 61 - 2a2, -bo)A{z) 

+ (bi - 02) [03^^ + a2Z^ + (60 - ai)z] , (3.29) 

and thus has two-step SI if and only if 60 = 2ai, as has been indicated by the second 
equality in (3.27). It is quite intriguing that a system has SI only when a parameter fulflUs 
a particular condition. This situation is reminiscent of conditional exact solvability [40]. 
Indeed, as was proved in the previous section, two-step SI always means solvability, and 
thus the present system is solvable only when bo = 2ai. That is, it is conditionally solvable. 
Hence, we shall call such SI conditional. 

The flrst and second intermediate potentials (3.21) and (3.22) in this case of conditional 
two-step SI read respectively as 

ii 15a3 5a2 , 4(6i)2 -h 30301 - (02)^ 

V [x; 61, bo) = ——z + -^ + —^7 ^- — ^z 

16 16 16(032;"^ -|- 022; -|- OiJ 



86160-0201 , 4(6o)2-(oi)2 



16(03^^ -|- 02Z -(- Oi) 16(03Z^ -(- O22; + Oi)z 



12 



(3.30) 

=z{x) 



and 



V'\x;bi,bo) = V'\x;bubo) - asz 
-2bo 



as&o ai{bi)'^ - 2a2bibo + 3as{bo 
bi bi{biz + bo) 



aiibif - aibibo + asibof 



b,{b,z + boY 



(3.31) 



z=z{x) 



Both of them do not have ordinary SI with V^ given in (3.28) irrespective of whether 
60 = 2ai or not. Hence, the conditional two-step SI of the system (3.28) is irreducible. 

The general solution to equation (2.12) with A{z) given in (3.27) is expressed in terms of 
elliptic functions when ai 7^ 0. The potentials are thus of deformed Lame or Heun type. We 
omit its involved general form here. When ai = 0, the system automatically has two-step SI 



since ftg =b 



Q — from the third equality in (3.27). In addition, the form of A{z) reduces 



to the one given in (3.36) with di = — a2, and the solution to (2.12) is given by (3.41). Sub- 
stituting (3.41) into (3.28), (3.30), and (3.31), we finally obtain the potentials in the x-space. 

When b^ 7^ — 2a2, we see from (3.24) and (3.25) that bi and bQ are functions of a^ 
(i = 0, . . . , 3) and thus are denoted by 



6- = 6f)-6f) = 2rfi(a). 



'(2) ,(0) 



^0 



^0 



2dt){a). 



(3.32) 



Substituting (3.32) into a general solution Qa^b^ = {b-^ +2a2)bf to the homogeneous equation 
(3.23), we obtain 



b'^' = '-^ibr'+d,)-do, 



3a 



(3.33) 



where and hereafter we omit the argument of di for the brevity. The equality (3.33) shows 
that bl and 6q are dependent parameters. Then, b\ and 6q must fulfill the same relation 
as theirs for the consistency. The latter requirement results in 



dn 



a2 + di 



du 



r(2) 



Q2 + di (2) 

"1 5 



r(0) 



3a3 3a3 

Substituting (3.34) back into (3.24) and (3.25), we have 



0-2 + di (0) 

~^ ^1 • 

6a^ 



ai 



«2 



{di 



ao 



(a2 + rfi)^(a2 - 2di] 



3a; ' """ 27(a3)2 

Under the latter parameter relations, the function A{z) is factorized as 



A{z) = a;\z + 



02 + di 



3a3 



z + 



02 — 2di 
3a; 



(3.34) 



(3.35) 



(3.36) 



and the pair of potentials (3.20) reads as 



V^{x-M):=V^{x-M.^^h, 



3a; 



ai = 



. (a2)^-(dl)^ 



ao- 



. (a2 + dl)^(a2-2dl) 



27(a3)^ 



3a3 a2±86i 3{di±2bi){3di±2bi) 
le"^ ^ 16 ^ 16(3a3Z + a2 - 2di) 



(3.37) 



z=z{x) 



13 



We now easily check that the latter pair of potentials satisfies 

V+{x- bi) = V-{x- bi + 2di) + 61 + rfi, 



(3.38) 



which is consistent with the formula for b\ in (3.32), and thus has two-step SI. The first and 
second intermediate potentials (3.21) and (3.22) in this two-step SI case read respectively 
as 



V'\x-M):=V''[x-M,''-^h, 



"1= 



(a2)^-(dl)2 „ (a2 + di)2(a2-2di) 



at) 



27(a3)^ 



ISag 5a2 ^{2bi + di){2bi - di) 
1^^ ^ Te" "^ 16(3a3^ + a2-2rfi) 



(3.39) 



z=z(x) 



and 



3(33 



"1- 3^; ' '^0- ,7f„„^2 



27(a3) 

as a2 + 16rfi 3(26i + rfi)(26i-cii) 
"16^ 48 ^ 16(3a3^ + a2-2rfi) 



(3.40) 



2=2(2:) 



Although both of them have the same functional dependence on the variable z as V^, neither 
has ordinary SI with V^. Hence, the system has irreducible two-step SI. The equation (2.12) 
with A{z) given in (3.36) is integrated as 



z + 



a-2 + di 
3a3 



di 1 2 / '^1 
— sech \ X 

a. V 2 



03 x^ 



for di 7^ 0, 
for di = 0, 



(3.41) 



Substituting (3.41) into (3.37), (3.39), and (3.40), we finally obtain the potentials in the 
x-space. 



Case 1-5: 04 7^ 

In the last, we shall consider the case of a fourth degree A[z). In this case, the pair of 
potentials (2.36) are given by 



v^ix-MM 



AA{z) 



3(03)^ , , \ 3 

30302 — 60401 =F 203O1 ± 804O0 z 



+ I 3(02)^ — I204O0 — 



303O1 3(03)^02 



2 

/o a 3(o3)2ai 
- 302O1 — 603O0 

\ 404 

3(oi)2 3(03)^00 



404 



=F4o26i ±60360 + {bif z 



4o4 



=F 80061 ± 2oi6o + 



404 

=F 60161 ± 40260 + 26160 ) z 

02 , 3(03)2 



I604 



±61 



z=z{x) 



(3.42) 
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The first intermediate potential (2.37) is calculated as 



V'\x;b,,bo) 



16A{z) 
-4(a2)' 

+ 4(6o)' 



80401 — 40302 + 



04 



z-^ + 4:{bi) + I604O0 + 203O1 



(03)^02 



04 



z'^ + I 86160 + 803O0 - 402O1 + 



(03)^01 



Oi 



(03)^00 



04 



10 2 I , "2 (03 

+ 2042; + 032; + — 



04 
^2 



I604 



(3.43) 



z=z{x) 



while the second intermediate potential (2.38) which exists only when 61 7^ is as 



V'\x; 61, 60) = V'\x; 61, 60) - 2o4^^ - 03^ - 



asbibp - 204(60)^ 



oi(6i)3 - 202(61)^60 + 30361(60)2 - 404(60) 



+ 2 



(61)2(61^ + 60) 
oo(6i)4 - oi(6i)36o + O2(6i)2(6o)2 - 0361(60)^ + 04(60)^ 



(6i)2(6i2 + 6o)2 



(3.44) 



z=z[x) 



Although V~^ and V~ have the same functional dependence on the variable z, it does not 
necessarily mean that there is a transformation of the parameters 61 and 60 which lead to 
non-refiective two-step SI. In fact, they satisfy the condition for two-step SI (2.3) if and only 
if 

/ O3 — 4o4 \ 



V 



as 

b^ + 4o2 -603 
6q + 601 b^ — 4o2 



5^0 



bt 
6+ 



0, 



(3.45) 



2oi/ 



where bf = b\ ± b\ . The trivial solution b^ =bQ = {i just produces reflective two-step SI 
(2.39). To have a non-trivial solution, the parameters must satisfy 

204(6^+402) = 3(03)2, (3.46) 

03(67 - 402) = -4o4(6(^ + 601), (3.47) 

03(6(^ - 2oi) = -3204O0. (3.48) 

We first consider the particular case when 03 = 0. In the latter case, the set of equations 
(3.46)-(3.48) reduces by the assumption 04 7^ to 

bt 
Hence, we have 

^(2) _ 1,(0) 



6^ + 4o2 = 60 + 601 = oo = 0. 



l(2) 



l(0) 



b^i' = ^i^"* ~ 4o2, 60^'' = —60^'' = — 3oi, A{z) = a^z + O22 + Oi^;. 
In this case, the pair of potentials (3.42) reads as 

1 



(3.49) 

(3.50) 



V^^(a;;6i,6o) 



-204(301 + 46o)z^ + (3o2 + 6i)(o2 + bi)z^ 

3(oi)2 



4:A{Z 

+ (30201 + 60161 ± 40260 + 26160)2 + 
+ 80061 ± 2oi6o + (6o)2 



y±6i. 



(3.51) 



15 



We can easily check that the latter pair satisfies 



V^{x- bi, bo)A{z) = V-{x; hi - 4a2, -bo)A{z) 

+ (61 - 2a2) [2042^ + 2a2z'^ + {bo - ai)z] 



(3.52) 



and thus has two-step SI if and only if bo = 3ai, as has been indicated by the second equality 
in (3.50). Hence, the present system has conditional two-step SI. 

The first and second intermediate potentials (3.43) and (3.44) in this case of conditional 
two-step SI read respectively as 



r\x-b,,bo) 



4:A{z) _ 
+ (bo)' 



204012^ -|- ((&! 



2aAZ^ + 



(02)^) z'^ + {2bibo - a2ai)z 
a-2 



(3.53) 



z=z{x) 



and 



V'\x- 61, 60) = V'\x- 61, 60) - 2a4z' + 



2a4(M' , ai(&i)=^- 202(61)^60-404(60)' 



{hY 



{b^Yib^z + bo 



26, 



01(61)^-02(61)260-04(60)' 



(6i)2(6i2 + 6o) 



(3.54) 



z=z{x) 



Both of them do not have ordinary SI with V^ given in (3.51) irrespective of whether 
60 = 3oi or not. Hence, the conditional two-step SI of the system (3.51) is irreducible. 

The general solution to equation (2.12) with A{z) given in (3.50) is expressed in terms of 
elliptic functions when oi 7^ 0. The potentials are thus of deformed Lame or Heun type. We 
omit its involved general form here. When oi = 0, the system automatically has two-step SI 
since 60 = 69 = from the second equality in (3.50). In addition, the form of A{z) reduces 
to the one given in (3.60) with 03 = 0, and the solution to (2.12) is given by (3.66). Sub- 
stituting (3.66) into (3.51), (3.53), and (3.54), we finally obtain the potentials in the x-space. 

When 03 7^ 0, we immediately obtain from (3.46) and (3.48) 



K 



-4o2 + 



3(03)= 



32O4O0 



2oi. 



(3.55) 



204 03 

Substituting them into (3.47), we have the following constraint 

64(04)^(40400 — 0301) = —(03)^ [I604O2 — 3(03)^] . (3.56) 

On the other hand, the homogeneous equations (3.45) under the conditions (3.46)-(3.48) 
lead to a single equation 036^ = 4o46||. Eliminating 6^ , 6q , and oq in it by using (3.55) 
and (3.56), we obtain 

3(03)3 



0361 



(0) 



Aa^bf 



I204O1 — 603O2 -|- 



204 



0. 



(3.57) 
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This relation means that bi and bo cannot be independent parameters. It is evident that it 
must hold not only between b[ and feg but also between bi and b^ for the consistency. 
From the latter requirement, we have another constraint: 



80401 — 40302 H = 0. 



04 



(3.58) 



Combining (3.56)-(3.58), we finally obtain 



bo = ^fci, b^? = bf^ - 4o2 + ^. 

404 /O4 



Oi 



0302 (03J 



Oo 



(03)^02 5(03)' 



(3.59) 



204 8(04)2' 16(04)2 256(04)3 ■ 

Under the latter parameter relations, the function A{z) is factorized as 



03 \ / 2 , «3 «2 5(03)2 



^^'^ = "^ < ^ + 4^ J 1,^ + 2^^ + ^ 16(04)2 



(3.60) 



and the pair of potentials (3.42) reads as 



V^ix-bi):=V^(x-bu^br') 



_ a3°2 ("3)^ 



(°3)^''2 S(a3)4 



2^4 8(a4)2' "^^ 16(04)2 256(04)8 



404 



16(04)2^2 + 8O4O32; + I6O4O2 — 5(03)^ 



3(«3)' , \ f 9(03)2 , \ a, 3(03)2 , 



3(03)^ 



. (3.61) 



z=z{x) 



We now easily check that the latter pair of potentials satisfies 



V+{x;b,) = V-{x;bi-Aa2 + ^^\ +2 



61 - 2o2 + 



4o4 



(3.62) 



which is consistent with the formula for b\ in (3.59), and thus has two-step SI. The first and 
second intermediate potentials (3.43) and (3.44) in this two-step SI case read respectively 
as 



il/__^ «3 



404 



0302 ("3)^ „ (°3)'^''2 5(g3)^ 
" 2a4 8(04)2 ' " 16(04)2 256(o4)3 



404 



16(04)2^2 -|- 8O4O32 -|- I6O4O2 — 5(03)2 

3(a3)2^'' 



X 



ib,r 



02 



804 



2 '^2 

+ 2042 + asz + — 



(«3)^ 



I604 



(3.63) 



z=z{x) 
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and 



V''{x;h):=V'Hx;b^,^h 



„,_ °3''2 (°3)^ „„— (<^3)^<^2 5(a3)"' 
"1- 2a4 ^(^^' "0 16(a4)2 256(a4)3 



4a/i 



16(04)^2;^ + 804032; + 160402 — 5(03)" 



X 



ihf - 02 



804 



2\ 2 



02 _ 3(03)^ 
2 I6O4 



(3.64) 



^=2(0;) 



Hence, the second intermediate potential V^"^ together with V^ constitutes a sequence of 
ordinary SI as 



V'\x; 61) = V- (x; h - 2o2 + -^) + 61 - 02 + 
V 404 / 

y+(x; h) = V'' (x; h - 2o2 + ^^) + 61 - 02 + 



3(03)2^ , . „ , 3(03)^ 

804 
3(03)2 



(3.65) 



804 



Therefore, the system has reducible two-step SI. The equation (2.12) with A{z) given in 
(3.60) is integrated as 



4yce- 



/2a4cx 



Z + 



Q3 
404 



a2\/2a4cx A 

1 



where 



'2aAX 



02 3(03)^ 



for c 7^ 0, 
for c = 0, 



(3.66) 



04 8(04)2 



(3.67) 



Substituting (3.66) into (3.61), (3.63), and (3.64), we finally obtain the potentials in the 
x-space. 

To summarize, we have found that a type A 2-fold SUSY system with an intermediate 
Hamiltonian resulting from a polynomial A{z) of at most fourth degree can possess non- 
reflective two-step SI, always when the degree is less than or equal to two, and under certain 
conditions when it is three or four. Another remarkable result is that all the type A 2-fold 
SUSY systems which have two-step SI as well in Cases 1-1 to 1-3 turn to have also ordinary 
SI, that is, they have reducible two-step SI, and they are all well-known SI potentials. On 
the other hand, all the other two-step SI potentials in Cases 1-4 and 1-5 are, to the best of 
our knowledge, new and except for the last model in Case 1-5 they do not possess ordinary 
SI, that is, they have irreducible two-step SI. In the next section, we shall examine more 
general cases where the function A{z) is not given by a polynomial. 

IV. NON-POLYNOMIAL A{z) 

To investigate the possibility of other irreducible two-step SI, let us coming back to the 
form of potentials (2.36). In general, the function A(z) can depend on the parameters bi. In 
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this case, the change of variables z = z{x) also depends on them via (2.12), and analysis of 
the two-step SI condition (2.3) with the pair of potentials (2.36) gets quite involved. In this 
work, we shall restrict our analysis to the case where A{z), and thus z = z{x) as well, does 
not depend on the parameters 6,. Under this restriction, we easily see from (2.36) that the 
two-step SI condition (2.3) holds if and only if the function A{z) satisfies 



4:A{Z) 



^(o)^2 



'(0) 



(0)^ 



'(0) 



4:A{Z) 



{h^z + Cr + 2((6r^ + O^'(^) - WA{z)) 



\b?z + b^^^Y - 2{{hfh + b^^^)A'{z) - 2b?A{z))] + m^^KbP) 



The latter condition is identical to the first-order linear differential equation 

{b+z + b+){b^z + 6o ) - 2(6|z + b^)A'{z) + 4RA{z) = 0, 



(4.1) 



(4.2) 



l(2) 



l(0) 



where bf = b\ ± b\ and R = R2 + bf. It is apparent that its most general solution A{z) 
depends on the parameters b\ . Hence, to obtain a solution without such dependence, we 
must impose the following additional condition: 



dA{z) 



db[ 



(0) 



0. 



(4.3) 



As we will show in what follows, the latter condition completely determines the admissible 
form of the function A{z) and the parameters bl and R2- 

To solve the differential equation (4.2), we first note that we can assume without any loss 
of generality that either b^ or 6^ is non-zero. Otherwise, we have 



'(2) 



'(0) 

h 5 



'(2) 



l(0) 
^0 ' 



i?2 = 0, 



(4.4) 



which exactly leads to the two-step SI of the type (2.39). Thus, bfz + b^ is not identically 
zero, and (4.2) is integrated as 



2A{z) = exp / dz 



2R 



dz {bl 2; -|- 6q ) exp 



dz' 



2R 



btz + btJ J ^'^'"''^ ' ^''^'^''K J ^'^btz' + b^ 
It requires separate treatments according to whether 6^ is zero or not. 



Case 2: 6+ 7^ 

Let us first study the case of non-zero bf . In this case, equation (4.5) reads as 



(4.5) 



2A{z) = [bfz + b+r / dz (b-,z + b-,){btz + 6+)- 



{btz + blY \ dz 



K 



{btz + btY-'^ + 



Kbo - b^b'o 



ibtz + b^y 



bt'' "" bt 

where fi = 2R/bi . Hence, we have the following three inequivalent cases. 
Case 2-1: /x^ 1,2 



(4.6) 
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In this case, we obtain from (4.6) 

A{z) = a2Z^ + aiz + oq + c{htz + h^^Y, (4.7) 

where c is an integral constant and 



(32 — 7777; 7; ^1 



2(2 -/i)' 2(l-/i)(2-/x)6+ 2(1 -/i)' 

Oo - 



(4.8) 



2(l-/i)(2-/i)(6+)2 2(l-/i)6r 
The condition (4.3) in this case is equivalent to the following set of equations: 
da2 dai da^ _^ dc dbf . dc db^ dfi 

We assume c 7^ to avoid duplication of Case 1-3. The last equality just means that /i does 
not depend on h\ . Substituting the expression of the parameters (4.8) into (4.9), we see 
that the following equations must hold: 

% = 0^ % = 0^ ^"^-^^^ = 0- (4-10) 

dhf^ dbf^ dhf^ dhf^ 

From the first and second equalities in (4.10), we immediately have 

6f ) = 6f) + 2rfi, hf^ = hf + 2rfo, (4.11) 

where di and d^ are constants which do not depend on h\ . 

Until now, we have not made any assumption on the relation between b\ and % , that 
is, they can be independent or dependent parameters. So, let us first assume that they are 
independent. In this case, we shall first show that b^ = 0. For this purpose, we note that 
so long as b^ = 26q + 2do j^ together with the assumption 6^ = 2b[ + 2di ^ 



dbf _ db+ _ db+ _ dbj 



dbf^ 4°) ^' dbf^ 4°) °- ^^-^^^ 



Then, the third equation in (4.10) results in b^ = b^ = 0, which is obviously contradictory. 
Hence, we must have fej" = 0, which together with (4.11) means 

6f ) = 6f) + 2du b^o^ = -bl^^ = do, (4.13) 

and the third condition in (4.10) is trivially satisfied. All the remaining equations in (4.9) 
to be satisfied now read as 

and their general solution is given by 

c = co/ibtr, (4.15) 
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where cq and /i are constants which do not depend on h\ . Substituting (4.13) back into 
(4.8), we have 

di = —{fi — 2)a2, do = —{fi — l)ai, ao = 0. (4.16) 

With the aid of the solutions (4.13), (4.15), and (4.16), we finally obtain the admissible form 
of the function A{z) and the parameters b] and R2 as 



A{z) = a2Z^ + aiz + cqz^, hf^ = bf^ - 2{fi - 2)02, 

bf^ = -4°) = -(/i - l)ai, R2 = {fi- 2)[6S°^ - (/x - 2)02]. 

The pair of potentials (2.36) in this case are given by 



(4.17) 






IQV^ix; 61, bo)A{z{x)) = 4[(a2)' + (61)'] z^ + 4(a2ai T 2ai6i ± 4a26o + 2&ifeo) 
+ 3(ai)2 ± 8ai6o + 4(6o)^ - 4 [(/i^ -Afi + 2)a2 T 2{fi - 2)bi\ cqz^ 
- 2^ [(2/i - 5)ai T 46o] cq^^"' - ^i{^i - ^){cofz^^-^\^^^^_^y (4.18) 

We can easily check that the latter pair satisfies 

V+{x- 61, b^)A{z) = V-{x- bi - 2(/i - 2)a2, -bo)A{z) 

+ [bi - (/i - 2)02] [(/i - 2)a2^^ + (&o - ai)z + (/i - 2)co2^] , (4.19) 

and thus has two-step SI if and only if 60 = (/^ — l)«'i, as has been indicated by the third 
equality in (4.17). Hence, the present system has conditional two-step SI. 

We note that 61 7^ in this case and the system always admits two intermediate Hamil- 
tonians. The first and second intermediate potentials (2.37) and (2.38) are respectively 
calculated as 

lW\x- 61, b^)A{z{x)) = A[{a2f + [bif] z^ + 4(a2ai + 2b^bo)z - {aif + 4(6o)' 

+ 4(/i2 - 2/i + 2)a2Co^^ + 2/i(2// - 3)aiCo^''^^ + /i(3// - 4)(co)';2'^"'|^^^(^), (4.20) 

and 

v'\x-MM) = v'\x-b,M + 
-2b 



\2i^.. ^.^ _^n\,^.v. i. \ , ai&i - 2a2&o - (/^ - 2)co6i2;^ ^ 



bxz + 60 
ai^i - 02^0 + c^biz^"^ 



{biz + bof 



(4.21) 



z=z{x) 



Both of them do not have ordinary SI with V"^ given in (4.18) irrespective of whether 
bo = {fi — l)ai or not. Hence, the conditional two-step SI of the system (4.18) is irreducible. 
We note that the present system (4.17)-(4.21) reduces to the one (3.27)-(3.31) in Case 1-4 
when we put /i = 3 with cq = a^ and to the one (3.50)-(3.54) in Case 1-5 when we put /i = 4 
with Co = 04. 

The equation (2.12) with A{z) given in (4.17) cannot be integrated analytically unless 
Oi = 0. When ai = 0, the system automatically has two-step SI since b^ = b^ = from 
the third equality in (4.17). In addition, the form of A{z) reduces to the one given in (4.24) 
with ai = 0, and the solution to (2.12) is given by (4.31). Substituting (4.31) into (4.18), 
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(4.20), and (4.21), we finally obtain the potentials in the x-space. 

In the next, let us consider the case when b{ and % are dependent. In this case, we 
can assume that % is a function of 6]^ . The set of conditions (4.9), and (4.10) as well, 
should be considered only for i = 1. The solution (4.11) to the first and second equations in 

(4.10) is still valid, but the one to the third one is b^ = zib^ where zi is a constant which 
does not depend onb[ . Substituting the latter into (4.11) and requiring that the functional 
dependence of 6q on b[ must be the same as that of 6q on bl , we obtain 

6(°) = zib^°\ do = zidi. (4.22) 

All the remaining equations in (4.9) are satisfied with the solution (4.15). Substituting 

(4.11) and (4.22) back into (4.8), we find that 02 7^ is necessary for the non-triviality of 
the system and obtain 

fef) = 5(0) _ 2(^ _ 2)a2, 2a2Zi = ai, 4a2ao = (ai)'. (4.23) 

With the aid of the solutions (4.11), (4.22), and (4.23), we finally obtain the admissible form 
of the function A(z) and the parameters b\ and R2 as 

A{z) = aJz + -^]\co(z + ^)\ 6f)=6f)-2(/i-2)a2, 

V 202/ V 2a2j (^4 24) 

"«' = i-/'' • "'"' = 1^"'°' ' «^ = (" - 2) ["'1°' - (" - 2)« J ■ 

The pair of potentials (2.36) in this case with the definition V{x; 61) := V{x; bi, 0161/(202)) 
are calculated as 

16V^{x■,b^)A{z{x))=4[{a2r + {b^f] (^ + ^) - 4 [(/.^ - 4/i + 2)a2 
T2(/i-2)6i]co('^ + ^) -/i(/i- 4) (00)^(^2; + ^^ 



2fi-2 

. (4.25) 

z=z{x) 



We can easily check that the latter pair satisfies 

V+{x; 61) = V-ix; 61 - 2(/i - 2)02) + {^i - 2)[6i - (^ - 2)02], (4.26) 

which is consistent with the formulas for b[ and R2 in (4.24), and thus has two-step SI. 
The first and second intermediate potentials (2.37) and (2.38), which always exist thanks to 
61 7^ 0, are respectively given by 

16V'\x-b,)Aiz)=4[ia2f + {bif] (^ + ^) + 4(/i2 - 2^i + 2)a2Co ("^ + ^ 

+ M3/.-4)(co)2(z + ^y , (4.27) 

and 

V'\x; 61) = V'\x- 61) - (/i - 2)co (z + ^\ . (4.28) 



2a2/ 
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Both of them in general do not have ordinary SI with V^ given in (4.25). Hence, the present 
system has irreducible two-step SI in general. It is worth noting, however, that the system 
(4.25)-(4.28) reduces to the reducible two-step SI system (3.61)-(3.64) in Case 1-5 when 
fj, = 4. Indeed, we can check that with the following substitution 

3(03)2 0302 3(03)3 

Co = 04, 02 = as -^, ai = — (4.29) 

804 2a4 16(04)^ 

the system (4.25)-(4.28) coincides with (3.61)-(3.64) with a^ replaced by Oj when /i = 4. 
On the other hand, when /i = 3, the system (4.25)-(4.28) reduces to the one (3.37)-(3.40) 
in Case 1-4 by the following substitution 

T 2di{a2 + di) , . 

Co = 03, ^2 = -di, ai = — , (4.30) 

keeping the irreducibility intact. Therefore, the reducibility of two-step SI in the system 
(3.61)-(3.64) in Case 1-5 is rather exceptional. 

The equation (2.12) with A{z) given in (4.24) is integrated as 

^ + ^V '' = ^tanh2(/i-2)./^x. (4.31) 

2a2/ 02 V 2 

Substituting (4.31) into (4.25), (4.27), and (4.28), we finally obtain the potentials in the 
x-space. 

Case 2-2: // = 1 

In this case, we obtain from (4.6) 

A{z) = a2Z^ + aiz + oq + cip^z + fe^) In \b'[z + 6^|, (4.32) 



where 



61 b,bt ., , , btbr, - b. 6(1 

J- ai = ^-f + cbt, flo = cb+ c = ^ " , ^ " 
2 ' 26+ ^ ° 2(6+ 



a2 = ^, ai = ^ + c6+, ao = c6+, c = "^ "°^ . "^ "° , (4.33) 



and c is an integral constant. The condition (4.3) in this case is equivalent to the following 
set of equations: 

da2 ^ d<M_ _db^ ^ dao_ _db^ ^ djcbt) ^ g(c6+) ^ 

dbf^ dbf^ dbf^ dbf^ dbf' dbf^ dbf^ ^^ 

We assume c 7^ to avoid duplication of Case 1-3. Substituting the expression of the 
parameters (4.33) into (4.34), we see that the condition (4.10) must hold also in the present 
case. So, let us first investigate the case when bl and 6q are independent. In this case, 
by following the same argument below (4.10), we arrive at the same solution (4.13). All the 
remaining equations in (4.34) to be satisfied now read as 



dc ^biC + do dc 



W^2 + -^7(0)-°' ^^-2^) 
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and their general solution is given by 

b+c = co-doln\bt\, (4.36) 

where Cq is a constant which does not depend on b\ . Substituting (4.13) and (4.36) back 
into (4.33), we have 

di = a2, ai = Co -(ioln|6||, Oq = 0, c = do/b^. (4.37) 

With the aid of the solutions (4.13) and (4.37), we finally obtain the admissible form of the 
function A{z) and the parameters b\ and R2 as 



(4.38) 



A{z) = a2z'^ + CqZ + dozln \z\, b\ = b\ + 2a2, 
6?) = -C) = c^o, R2 = -^r - a2. 
The pair of potentials (2.36) in this case are given by 

16\/^(x; 61, bo)A{z{x)) = 4 [(02)^ + {bif] z^ + 4(a2Co + 202^0 T 2co6i ± 2dobi 
± 40260 + 26i6o)^ + 3(co)2 + 2corfo + 3(^0)^ ± 8(co + do)bo + 4:{bof 
+ 4do(a2 T 2bi)z\n \z\ + 2cio(3co + do± 46o) In \z\ + 3(cio)^(ln \z\f\^^^^^y (4.39) 

We can easily check that the latter pair satisfies 

V+{x; 61, bo)A{z) = V-{x; b, + 2a2, -bo)A{z) 

- (61 + 02) [022^ + (cq - do - bo)z + doz In |2;|] , (4.40) 

and thus has two-step SI if and only if bo = —d^, as has been indicated by the third equality 
in (4.38). Hence, the present system has conditional two-step SI. 

We note that 61 7^ in this case and the system always admits two intermediate Hamil- 
tonians. The first and second intermediate potentials (2.37) and (2.38) are respectively 
calculated as 

16r^(x; 61, bo)A{z{x)) = 4[(a2)' + (&i)'] z^ + 4(a2Co + 2bibo)z - (cq - dof 

+ A{bof + 4a2dozln \z\ - 2do{co - do) In [^l - {dof{\n \z\f\^^^^^y (4.41) 

and 

12/ 1. I. \ T/ii/ I. I. \ , {co - do)bi - 2a2bo + dobiln\z\ 



V'\x;bubo) = V''{x;bubo) + 
-26, 



61-2; + bo 
cobi — 02^0 + dobi In \z\ 



ipiz + 60) 



(4.42) 



z=z{x) 



Both of them do not have ordinary SI with V^ given in (4.39) irrespective of whether 
60 = — c^o or not. Hence, the conditional two-step SI of the system (4.39) is irreducible. Un- 
fortunately, the equation (2.12) with A{z) given in (4.38) cannot be integrated analytically 
in general for rfg 7^ 0. 
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In the next, let us consider the case when b\ and % are dependent. In this case, 
the general solutions to (4.10) are given by (4.11) with (4.22). However, with the latter 
solutions c will be zero by the last equality in (4.33), and we have a quadratic A{z) as a re- 
sult. Hence, the system with the dependent bi and Bq just reduces to a system in Case 1-3. 

Case 2-3: /i = 2 

In this case, we obtain from (4.6) 

A{z) = a2Z^ + aiz + ao + cib^z + b^fhi \b'[z -h &(j"|, (4.43) 

where 

bib, - b-,bt 



0-2 = c{bj) , ai = 2cb{b^ - 



(4.44) 



and c is an integral constant. The condition (4.3) in this case is equivalent to the following 
set of equations: 






dc _db'{ _^ dc _(96q 

W^ W^ 'W^ d^^ 



bt-M + 2c:Jm = K-m^ + 2C3 = 0. (4.45) 



We assume c 7^ to avoid duplication of Case 1-3. Substituting the expression of the 
parameters (4.44) into (4.45), we see that the condition (4.10) must hold also in the present 
case. So, let us first investigate the case when bi and fog ^i^e independent. In this case, 
by following the same argument below (4.10), we arrive at the same solution (4.13). All the 
remaining equations in (4.34) to be satisfied now read as 

db?^ ibtr 'db^ ' ^^ ^ 

and their general solution is given by 

(6+)2c = Co-rfiln|6+|, (4.47) 

where Cq is a constant which does not depend on b\ . Substituting (4.13) and (4.47) back 
into (4.44), we have 

a2 = cq — di\n\bf\, ai = —do, ao = 0, c = di/(bfY. (4.48) 

With the aid of the solutions (4.13) and (4.48), we finally obtain the admissible form of the 
function A{z) and the parameters bl and R2 as 

A(z) = cqz'^ + aiz + diz^ In \zl bf^ = 6S°^ + 2di, 
% = -% = -a-i, i?2 = 0. 
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The pair of potentials (2.36) in this case are given by 

16V^{x; bi)A{z{x)) = [4(co)^ + 3(t/i)2 ± Sdih + 4(6i)2] z^ + 2[ai(2co - 3di) 
T 4:aibi ± 4(2co + di)bo + Abibo\z + 8(01)^ ± Sai^o + 4(6o)^ 
+ Scodiz'^ In \z\ + 4di(ai ± Abo)z In |2| + 4(cii) V(ln \z\f\^^^^^y (4.50) 

We can easily check that the latter pair satisfies 

V+{x; 61, bo)A{z) = V-{x; bi + 2d,, -bo)A{z) + {bo - ai)(&i + di)z, (4.51) 

and thus has two-step SI if and only if bo = ai, as has been indicated by the third equality 
in (4.49). Hence, the present system has conditional two-step SI. 

We note that 61 7^ in this case and the system always admits two intermediate Hamil- 
tonians. The first and second intermediate potentials (2.37) and (2.38) are respectively 
calculated as 



16V"'^(x; bi, bQ)A{z{x)) = [4(co)^ + Scod, - {d 
+ 2(2aiCo + baidi + 4:bibo)z - {a{f + 
+ 8(ii(co + di)z^ In \z\ + Aaidizln. \z\ + 



i?+Kb,fy 


^ihof 


-A{d{)''z'{\xi zf 



--z(xy 



(4.52) 



and 



y'\x-MM 



V'\x-MX)-do 

a\b\ — Co^o ~ c^i^o In k 



a\b\ — (2co — d\)bo — 2d]bo In \A 



2bn 



{biZ + bo 



biz + bo 



z=z{x) 



(4.53) 



Both of them do not have ordinary SI with V^ given in (4.50) irrespective of whether bo = cti 
or not. Hence, the conditional two-step SI of the system (4.50) is irreducible. 

The equation (2.12) with A{z) given in (4.49) cannot be integrated analytically unless 
Oi = 0. When oi = 0, the system automatically has two-step SI since 6q = 6q = from 
the third equality in (4.49). In addition, the form of A{z) reduces to the one given in (4.55) 
with fli = 0, and the solution to (2.12) is given by (4.59). Substituting (4.59) into (4.50), 
(4.52), and (4.53), we finally obtain the potentials in the x-space. 

In the next, let us consider the case when b\ and 6q are dependent. In this case, the 
general solutions to (4.10) are given by (4.11) with (4.22). All the remaining equations in 
(4.9) are satisfied with the solution (4.47). Substituting (4.11), (4.22), and (4.47) back into 
(4.44), we find that 02 7^ is necessary for the non-triviality of the system and obtain 



^2 = Co — di In \bf 



2a2Zi = fli. 



4a2ao 



ai 



(4.54) 



With the aid of the solutions (4.11), (4.22), and (4.54), we finally obtain the admissible form 
of the function A{z) and the parameters b\ and R2 as 



A{z) 



'(2) 



2; + -^ ) ( Co + di In 



z + 



2ao 



2a2 ' 



'(0) 



^bf\ 
2a2 ' 



R2 



(2) 



bf^ + 2rfi, 



(4.55) 
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The pair of potentials (2.36) in this case with the definition V{x; 61) := V{x; &i, ai&i/(2a2)) 
are calculated as 



16V^{x;bi)A{z{x)) = (z + 



2ao 



+ Scodi In 



4(co)' + (c/i±26i)(3rfi±26i; 



z + 



2^ 



4(rfi)^ In 



z + 



2^ 



(4.56) 



2=2(2:) 



We can easily check that the latter pair satisfies 

V+{x-h) = V^{x-bi + 2di), 



(4.57) 



which is consistent with the formulas for b[ and R2 in (4.55), and thus has two-step SI. 
The first intermediate potential (2.37) is given by 



16V'\x; b^)A{z)= (^ + ^) 



^cof + 8cod,-{d,Y + ^b^ 



+ 8(co + di)di In 



z + 



2ao 



+ 4(rfi)Mln 



z + 



2a2 



(4.58) 



z=z{x) 



Intriguingly, the second intermediate potential (2.38) which always exist thanks to 61 7^ 
differs from the first one only by a constant, V^(x; 61) = V^^{x; bi) — di. Both of them do 
not have ordinary SI with V^ given in (4.56). Hence, the present system has irreducible 
two-step SI. 

The equation (2.12) with A{z) given in (4.55) is integrated as 



In 



z + 



2^ 



di 2 _ '^0 
~2^ "(fi 



(4.59) 



Substituting (4.59) into (4.56) and (4.58), we finally obtain the potentials in the x-space. 



Case 3: 6+ = 

In this case, we obtain from (4.6) 



2A{z) = e"' f dz {-2bfz + b^)e 



(4.60) 



where u = 2i?/feJ". We can assume without any loss of generality that z/ 7^ 0. Indeed, if 
i^ = 0, we have 



(0) 2 



2A{z) = -br z^ + bo z + 2c. 



(4.61) 



where c is an integral constant, and thus it just reproduces the case of quadratic A{z). When 
z/ 7^ 0, we obtain 



A(z) = aiz + ao + ce'' 



(4.62) 
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where c is an integral constant and 

bf^ _2hf-h,v 

The condition (4.3) in this case is equivalent to the following set of equations: 

(9a 1 da^ dc du 



ai = ^, ao = ' ,,,, ° ■ (4.63) 



dbf^^ dh?^ dh?^ dh?^ 



0. (4.64) 



The last two equalities just mean that both c and u do not depend on 6^ . We assume c 7^ 
to avoid duphcation of Case 1-2. Substituting the expression of the parameters (4.63) into 
the first two equalities in (4.64), we have 



dhf^ _ dh_ 



0, -^ = 0. (4.65) 



The first equation in (4.65) is inconsistent if hi is an independent variable. Thus, the system 
has only one independent variable h^ , and h\ is a. function of it. Then, we immediately 
have h\ = di where di is a constant which does not depend on 6q . The same equality 
must hold also for b[ for the consistency, bl = di. On the other hand, bf = b\ +bi = 
in the present case and thus b[ = b[ = is the only permissible solution. The general 
solution to the second in (4.65) is obviously fep = ^0 + ^c^o where do does not depend on 
bg . Substituting these solutions back into (4.63), we obtain 

ai = 0, do = —aou. (4.66) 

With the aid of the latter solutions, we finally obtain the admissible form of the function 
A{z) and the parameters b\ and R2 as 

A{z) = ao + ce"^ bf = bf - 2aoz/, R2 = ipf - aov)v. (4.67) 

The pair of potentials (2.36) in this case are given by 

\W^{x- bo)A{z{x)) = -cV^e^- - Aci^iaou t 26o) e^^ + ^bof\^^^^^y (4.68) 

We can easily check that they satisfy the two-step SI condition (2.3) with the parameters 
bo and R2 given in (4.67). In this case, 61 = and the system admits only one intermediate 
Hamiltonian. The intermediate potential (2.37) is calculated as 

16r^(x; bo)A{zix)) = ^c^u^e^"' + Aaocu^e"' + 4(6o)'. (4.69) 

Although it has the same functional dependence on the variable z as V^ but is not ordinary 
SI. Hence, the system is irreducible two-step SI. 

The equation (2.12) with A{z) given in (4.67) is integrated as 



ao I tanh J —vx — 1 ) for Oq 7^ 0, 



ce^^={ \. V ^ J (4.70) 



y'ir^l 



for ao = 0. 



Substituting (4.70) into (4.68) and (4.69), we finally obtain the potentials in the x-space. 
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Case 1-1 


Case 


1-2 


Case 


1-3 


Case 


1-4 


Case 


1-5 


Case 


2-1 


Case 


2-2 


Case 


2-3 



Reducible 
Reducible 
Reducible 

Irreducible, Conditional (b^ = —2a2) 
Irreducible (6j~ 7^ —202) 
Irreducible, Conditional (03 = 0) 
Reducible (03 7^ 0) 

Irreducible, Conditional (2 independent parameters) 
Irreducible (1 independent parameter, ^ 7^ 4) 
Irreducible, Conditional 

Irreducible, Conditional (2 independent parameters) 
Irreducible (1 independent parameter) 
Case 3 Irreducible 

TABLE I: List of the obtained two-step SI potentials. Cases 1-4 and 1-5 can be included in Case 
2-1 (/x = 3 and // = 4, respectively). 



V. DISCUSSION AND SUMMARY 

In this paper, we have studied and constructed two-step SI potentials by utilizing the 
framework of A/"- fold SUSY. Recognizing the crucial fact that two-step SI always means type 
A 2-fold SUSY with an intermediate Hamiltonian, we have successfully revealed its general 
aspects and made the systematic construction of such systems. The essential point of our 
analysis resides in the fact that type A A/'-fold SUSY systems possess the simple general 
form (2.36) in terms of the variable z. It enables us to make the model-independent analysis 
without recourse to any specific assumption or ansatz. Furthermore, we have obtained even 
such two-step SI systems that admit an analytical expression only in terms of z, namely, in 
Case 2-1 with ai 7^ 0, Case 2-2, and Case 2-3 with ai 7^ 0. In Table I, we summarize the 
inequivalent two-step SI systems and their properties obtained in this paper. 

It is remarkable that all the obtained two-step SI as well as ordinary SI potentials are 
of translational classes, that is, parameter relations are characterized by constant shifts 
^i ^ ^i + 2'ii- In the traditional approaches, one assumes the latter relation from the 
beginning and then tries to construct an SI potential which meets it. A significant point of 
our analysis is that it is a consequence of the differential equation dh~ /dh] = originated 
from the requirement (4.3) that A{z) should not depend on the parameters hi ■ It in 
particular means that another class of SI potentials which is different from translational 
ones could be obtainable only from such a A{z) that depends explicitly on h\ . Then, to 
obtain a scaling SI potential which has parameter relations h\ = qib\ , for instance, SI 
conditions are expected to lead to a differential equation like 



A^(ln6f)-ln6f) = 0. 
dby 

In this way, we would be able to make a systematic analysis also on non-translational SI 
potentials. 
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In [14], the concept of multi-step SI was also introduced as a natural extension of two-step 
SI. Following a similar argument in Section 11, we easily conclude that multi-step SI is a 
special case of A/'-fold SUSY with intermediate Hamiltonians at every intermediate positions 
for M > 2. Until now on, there has been only one systematic study on such systems in the 
case of type A 3-fold SUSY [38]. It was shown in the latter reference that every type A 
3-fold SUSY systems with intermediate Hamiltonians at every intermediate positions, which 
were dubbed Class (1,1), were ordinary SI as well. That is, any three-step SI in type A 
3-fold SUSY is reducible. Hence, irreducible three-step SI can be realized, if it exists, only 
in other types of 3-fold SUSY. 

Another interesting aspect of our present results which can be immediately read from 
Table 1 is that all the conditional two-step SI are irreducible. To the best of our knowledge, 
there have been no decisive mathematical understanding of conditional exact solvability. 
The result thus indicates the possibility that the latter concept might be well characterized 
in a context of irreducible two- and multi-step SI. 

Although we have restricted our investigation to ordinary scalar Schrodinger operators, 
the concepts of ordinary, two-step, and multi-step SI are applicable to much wider sys- 
tems. In fact, TV- fold SUSY was successfully formulated also for Schrodinger operators 
with position-dependent mass [35], matrix ones [41], and ones with reflection operators [42]. 
Hence, we would be able to make systematic studies on various SI in these systems with the 
framework of A/'-fold SUSY, as have been done in this work. 
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